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We study the dispersion mechanism of Lieb mode excitations of both single and multi component 
ultra-cold atomic Bose gas, subject to a harmonic conhnement through chirp management. It is 
shown that in some parameter domain, the hole-like excitations lead to the soliton’s negative mass 
regime, arising due to the coupling between chirp momentum and Kohn mode. In low momenta 
region the trap considerably affects the dispersion of the grey soliton, which opens a new window to 
observe Lieb-mode excitations. Further, we extend our analysis to binary condensate, which yields 
usual shape compatible grey-bright soliton pairs. The inter-species interaction induces a shift in 
the Lieb-mode excitations, where the pair can form a bound state. We emphasize that the present 
model provides an opportunity to study such excitations in the low momenta regime, as well as the 
formation of bound state in binary condensate. 


I. INTRODUCTION 

Since the first successful realization of the atomic Bose- 
Einstein condensate (BEC), the experiments, as well as 
theoretical works on this unique quantum state has ad¬ 
vanced greatly by exploring different intriguing domains 
and expanding new horizons. By now we have obtained 
great controllability and tunability in the magneto-optic 
setup to cool, preserve and tune the atomic gases. This 
unique controllability has enabled multi facet research in 
ultra-cold atomic gases. Among them the quasi one di¬ 
mensional (cigar shaped) condensate and associated soli¬ 
ton solutions of the mean-field Gross-Pitaevskii equation, 
are quite well studied till date 

Nevertheless, the scope to study these non-linear sys¬ 
tems are never ending HHin] which has been substan¬ 
tiated through several investigations during the last few 
years. Among them, we note the study of a trapped BEC, 
under the infiuence of an oscillating Gaussian potential, 
where both vortex pairs and solitons can be created by 
suitable amplitude modulations m- In another work, 
splitting of the ground state of an attractively interact¬ 
ing BEC into two bright solitons with controlled relative 
phase and velocity has been explored [12]. The nonlin¬ 
earity of Gross-Pitaevskii (GP) equation, describing BEC 
at the mean-field level, leads to these excitations HMS]. 

One may also recall that in 1963, Lieb discovered a 
collective excitation of the condensed bosons in a sec¬ 
ond quantized formulation, which exhibited a periodic 
dispersion, very different from the well-known Bogoli- 
ubov mode m- The same was later identified with a 
complex grey soliton at the mean field level m, which 
mimics a hole-like excitation of the condensate. Labora¬ 
tory confirmation of this mode had been difficult due to 
the fact that, in the experimentally accessible low mo¬ 
menta regime, Lieb mode dispersion lies below that of 
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the Bogoliubov excitations na na [hhui- The insta¬ 
bility of BEC for higher momentum values, where these 
two modes differ significantly, excluded this domain from 
observation, for quite some time. 

Recent observation of oscillations between stable local¬ 
ized grey solitons in a cigar shaped BEC and vortex rings 
has led to significant interest in the Lieb mode, having 
complex envelope profile m- These solitons have been 
produced through collision of two BECs Eain, when 
the collisional energy is reduced to a level, where it is 
comparable with the interaction energy. The wavelength 
of the laser, responsible for density modulation, is larger 
than the healing length. In this domain, nonlinearity 
plays a stronger role than dispersion and the interference 
effect leads to the formation of an array of grey solitons. 
At the lowest collisional energy, the interference pattern 
produces a single pair of grey solitons. These coherent 
structures may have relevance for atom interferometry 
[S] I 23 H 25 ] . The presence of a harmonic trap and the time 
dependent coupling parameters, can non-trivially affect 
the soliton dynamics. In single component BEC, soliton 
and soliton trains in a trap with time varying scattering 
length and loss/gain have been analyzed [26]. A num¬ 
ber of recent studies have explored the dark and bright 
solitons in the cigar-shaped BEC, revealing the nontriv¬ 
ial coupling between trap and soliton dynamics [26H^ . 
Hence, it is natural to inquire about the nature of the 
coupling between the trap and soliton geometry. The 
key question that arises here is about the possibility of 
accessing the regime of Lieb-mode excitations, associated 
with a quantum grey soliton through the coupling be¬ 
tween trap geometry and their excitations. 

The aim of this letter is to address the above question 
and how the confining potential under certain circum¬ 
stances affects the hole-like Lieb mode excitation. To 
start with, we show exact grey soliton solutions for one 
and two component BEC subject to a harmonic confine¬ 
ment. It is interesting to observe the dramatic collapse 
and revival of the atomic condensate by appropriate tai- 
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loring of the loss/gain term. Further investigations re¬ 
veal that these chirped grey solitons can even be ac¬ 
celerated, compressed or brought to rest, when the ef¬ 
fects of the harmonic trap on the width and amplitude 
of the solutions are explored. We then study the soli- 
tonic Lieb-mode dispersion and compare it with the usual 
sound (Bogoliubov) mode dispersion. It is interesting 
to observe that the temporal modulation of the chirping 
is found to be important to resolve the Lieb mode ex¬ 
perimentally. The time modulated chirping enables us 
to find a small window where the hole-like excitation 
mode dominates over the sound mode before destabi¬ 
lizing the system. When the time dependency is over¬ 
looked, we extract the usual dispersion profile as shown 
in m- One must note that, our formalism also elab¬ 
orates on the control mechanism of chirping by means 
of external magnetic field. Thus the current formulation 
paves the way for their coherent control and manipu¬ 
lation at low energy, through the temporal modulation 
of the scattering length, as well as the trap [26l I321I38] . 
The two component BECs (TBECs) show a rich variety 
of solutions (dark-dark, dark-bright and bright-bright), 
though in this paper, we restrict ourselves to the dark- 
bright pair. The effects of the inter-species interaction 
and harmonic trap on the dispersion of these hole-like 
excitations are investigated. In case of attractive inter¬ 
actions, the energy shows an upward shift with respect 
to the single component case. As we move from attrac¬ 
tive regime towards repulsive, their energy decreases and 
exactly matches with the single component case, when 
inter-species interaction becomes zero. In the repulsive 
domain, the energy goes down further, entering into the 
negative regime, which leads to the formation of a bound 
state for these hole-like excitations. We argued that this 
bound state formation can manifest in the low energy 
mean-field theory. 

The letter is organized in the following way, in Sec. 
2 we present the exact analytical solutions for both one 
and two component case. These solutions lead to hole- 
like excitations or quantum grey soliton. We further show 
the dynamics in the presence of both regular and expul¬ 
sive harmonic trap in Sec. 3, where these excitations 
undergo nonlinear compression in some appropriate pa¬ 
rameter regime. We investigate their dispersion relation 
in the later section (Sec. 4) and analyze the effects of 
the trap and interactions on the hole-like excitation as¬ 
sociated with the Lieb mode. We draw our conclusion in 
Sec. 4 and discuss about the possible future goals. 


II. QUANTUM GREY SOLITON 
Single Component 

The dynamics of BEG is well described at the mean- 
field level by the following 3D GP equation, with the 


order parameter T: 

= - — + + ( 1 ) 

Here, V{r,t) = V{x,y) + V{z,t) and U = A.'nW'aglm 
is the strength of the atom-atom interaction, with 
and m being the scattering length and the mass of 
the atoms, respectively. At sufficiently low tempera¬ 
tures, the BEG confined in a strong transverse trap, 
V{x,y) = with uj±_ being the transverse 

frequency, can be made effectively one dimensional. 

The appropriately scaled quasi-ID GP equation, in di¬ 
mensionless units, can be written in the form [26l[35] . 

1 2 / 2 / I / 

^~dt ^ '<P+2^{t)z 

Here, the strength of the atom-atom interaction and 
the spring constant are 7 (t) = 2as{t)/aB and M(t) = 
respectively, with being the Bohr radius. 
Eor the sake of generality, we have kept M{t) time de¬ 
pendent. Regular and expulsive oscillator potentials cor¬ 
respond to M > 0 and M < 0, respectively. 

In order to solve the above nonlinear equation, a self¬ 
similar method is adopted here, which is straightforward 
and as a result of it, one can clearly see the dependence of 
the trapping potential and the time dependent nonlinear¬ 
ities. This method has already been elaborated in detail 
by Atre et. al [26]. With this approach, we assume the 
following ansatz solution: 

^p{z, t) = , (2) 

where, we described the solutions in terms of the cen¬ 
ter of mass co-ordinate = A{t){z — l{t)) and 

G{t) = g{t')dt' represents the phenomenological 
loss/gain term. The kinematic phase has a quadratic 
form: ^{z,t) = a{t) + b{t)z — ^c{t)z‘^^ exhibiting chirp¬ 
ing, where a{t) = uq — and h{t) = A{t). 

Here p = —2Kao = g X; fi = 2u(t)/A^(t) is the chem¬ 
ical potential in the moving co-ordinate frame and A, 
respectively, are the chemical potential and a constant 
parameter controlling the energy of the excitation. The 
coefficient of chirping c(t) can be determined from the 
following Riccati type equation, — c^{t) = M{t). It 
is worth pointing out that the above Riccati equation can 
be expressed as a Schrodinger eigenvalue problem by in¬ 
corporating the appropriate transformation as shown in 

m- 

The width of the soliton profile is given by, A{t) = 
AqcJ^o At )dt ^ worth noting that center of mass mo¬ 
tion (GOM) has a direct connection with the Kohn theo¬ 
rem [39] . Eor better understanding, we separate the gov¬ 
erning equation of GOM: -h c{t)l{t) = A{t)u. The 

Kohn theorem states that the center of mass motion os¬ 
cillates exactly with the trap frequency, undergoing sinu¬ 
soidal oscillations. When the trap frequency varies with 
time, the GOM motion ceases to decouple from the trap. 
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FIG. 1. Collapse and revival of the atomic condensate, 
(a) depicts the collapse and revival with increasing amplitude 
by appropriate tailoring the loss/gain = 0.1). (b) repre¬ 
sents the same collapse and revival of the atomic condensate 
with decreasing amplitude for a different value of loss/gain 
(/3 = — 0.1). In both the cases, the grey soliton shows small 
oscillation in space due to the presence of the harmonic con¬ 
finement. The other parameter values are as follows: Aq = 1, 
ao = 1.5, lo = 0.7 and A = 0.9, 70 = 0.5 and M{t) = 1. 


leading to a nontrivial behavior as discussed in detail 
by the authors in uni- In order to gain complete con¬ 
trol of the spatio-temporal dynamics of the condensates 
and these hole-like excitations, we separate the resulting 
equation of the COM of these excitations, 

^+M(i)^(i)=0. (3) 

From current conservation, amounting to solving the 
imaginary part of the GP equation, one obtains, = 

u{l — where, ao is the Thomas-Fermi background 

density of the atoms in the center of mass frame 
In this frame, the density equation can be cast in the 
convenient form [15], 


\ ) 


{ncF — V?) 


(g - gp)^ 

2g ^ 


( 4 ) 


where, 7 (t) = 7 oe with The solution 

of Eq. (pj) takes the self-similar form. 


r(^) = go — go cos^ 0 sech^ 


cos 0 




( 5 ) 


where the Mach angle is given by, 0 = sin ^ ^ = 

sin“^ Here, V = A(t)u is the velocity of the soli- 
tons, which is bounded by the sound velocity Cg = 
A{t)yJn(7Q = A{t)Cs] Cg is the sound velocity in the ab¬ 
sence of the harmonic trap. By appropriate tailoring of 
the loss/gain {g{t) = —Pt) leads to an interesting phe¬ 
nomena: collapse and revival of the atomic condensate. 
In Fig.0, the collapse and revival with decreasing &) 
and increasing &) amplitudes of the quantum grey soli¬ 
ton are shown for positive and negative loss/gain, P = 0.1 
and /3 = —0.1, respectively. The collapse and revival of 
this hole-like excitations are connected with the inter¬ 
play between trap geometry and time modulated inter¬ 
actions, which can be tailored through the loss/gain term 


as well. In the limiting case, i.e., when the harmonic trap 
is switched off, all the expressions match with the known 
results nsms]. 


Binary condensate 

Binary mixtures of BECs are available in experimental 
studies. Most of the binary mixtures contain two differ¬ 
ent hyperfine states of the same atomic species, such as 
[41], and ‘^^Na [42]. A BEC has also been cre¬ 
ated in a hetero-nuclear mixture of and [43] . 

The two component BECs (TBECs), present a novel and 
fundamentally different scenario for their ground state 
structure, as compared to the single component case. 
Theoretically, TBECs have been widely studied and var¬ 
ious structures have been identified, such as dark-dark 
[44] and dark-bright solitons [45]. Very recently, Ham- 
ner et.al.^ presented the first experimental observation of 
dark-bright soliton trains in superfiuid-superfluid coun¬ 
terflow, and showed that this induced a modulational 
instability in this miscible system [46]. In another paper, 
Shimodaira et. a/., considered two component BEC ro¬ 
tating in a toroidal trap and showed the transition from 
immiscible to miscible condensate m- 

In this analysis, we consider the binary condensate, 
placed in two different hyperfine states of ^^Rb atoms. 
The TBEC are described by the coupled GP equations, 
which can easily be put into dimensionless form, similar 
to the single component: 

k=a,b 

- 

where, ^j{zR) with j = a^b refers the wave function 
of the first and second component, respectively. The 
intra-species interactions between the atoms are given 
by, = 2^^^^, which in fiber optics language is 

known as self-phase modulation, while the inter-species 
interactions are denoted as, ^ab{t) = 2^^^, known as 
cross-phase modulation. Without any loss of generality, 
we here assumed jab{t) = ^ba{t) and from now onwards 
for notational convenience, we denote Ajj{^) — 
yj{t) = iJijA^{t) the chemical potential of component. 
The BEC wave functions are taken in the form, 

(7) 

where, the parameters i{zR)^ A{t) and c(t) are same 
as in the previous case, with l{t) = j^v(t')dt'. In 

this case, Gj(t) = gj{t')dt' is the phenomenological 
loss/gain term for first {j = a) and second component 
{j = b). Eor the sake of simplicity, we have considered 
here: a{t) = b{t) = 0. Therefore, from current conser¬ 
vation, we obtain the following two hydrodynamic equa- 


















4 



FIG. 2. The figures show the variation of the soliton profile 
for ^ = 7r/16. Fig. (a) and (b), respectively show the density 
profile of the hole-like excitations for = 0 and (5a — 0.05. 
The bright soliton prohle is shown in (c) and (d) for /3b = 0 
and /3b = 0.2, respectively. We consider here, = 4, bCb = 0.5 
and Ka — 2.4. The other parameter values are same as in Fig. 

0- 


tions, 


dXa 

di 


0-0 dxb 

uil -) and = u. 

o-a 


( 8 ) 


The phase of the second component is taken to be in¬ 
dependent of its density, implying that the correspond¬ 
ing superfluid velocity is constant. The center of mass 
motion of the soliton in this case, turns out to be: 
l{t) = loe~ ^0 , which can be controlled by modulat¬ 

ing the trapping frequency. The real part of the coupled 
equations leads to the following soliton solutions of the 
first and second condensates respectively, 


cTog) = CTo(l - cos^6»sech^(^^^^(2:,t)^, (9) 

<^biC) = ( 2 :, i)), (10) 

where, Ni, = f dz\ 2 pi){z,t)\‘^ is the rescaled number of 
atoms in state b. The interaction parameters are scaled as 
: K,j = with j = a, 6, ab. Akin to the single compo¬ 

nent case, this shows that the frequency of the harmonic 
oscillator is related to the interaction parameter through 
the chirping term. The chirped phase in this case plays 
a very crucial role in determining the various parameters 
associated with soliton s. The healing length is fou nd to 
be of the form: C = \^{Kab - K6)/[cro(K^j - «;»«&)], and 
9 = sin“^ —. The soliton’s velocity u achieves a maxi- 
mum value = ^, when 0 = ^. This corresponds to 


a constant solution a a = ctq. The 6> = 0 limit leads 
to a pure dark soliton solution aa = ao tanh[^( 2 (, t)/C], 
which is static, since rt = 0. The chemical potential 
of the first component is similar to the single compo¬ 
nent case, whereas, for second component it strongly de¬ 
pends on the inter-species interaction and healing length : 
i^b = i^abO'o — \v? — *^ 2 ^ 2 ^ • Similar to the single component 
case, we assume here the loss/gain term as gj{t) = Pjt. 

Fig. § depicts the corresponding density profile of 
the grey and bright solitons. Fig. &) shows the propa¬ 
gation of a grey soliton, which oscillates with time due to 
the presence of the harmonic trap, without the presence 
of loss/gain (^5 = 0). The phase associated with this exci¬ 
tation changes at each oscillating point. Fig. if) refers 
to the density distribution of the hole-like excitation with 
loss/gain (/3 = 0.05). The amplitude of the atomic distri¬ 
bution diminishes with time. The density profile of the 
bright soliton, characterized by a localized maximum is 
shown in Fig. ([^). It is clear from the figure that bright 
soliton also oscillates at the same frequency as in the case 
of grey soliton. This frequency of oscillation is directly 
related to the frequency of the harmonic trap. Fig. ®) 
depicts the same in the presence of loss/gain (/3 = 0.2). 
The amplitude of the bright soliton decreases with time 
as well, though the frequency of the oscillations remains 
same. One can vary the soliton profile by changing the 
nonlinearities and the trap frequency. In the real experi¬ 
ment, the nonlinearities and the confining harmonic trap 
can be changed by tuning the magnetic field induced Fes- 
hbach resonances and optical trap respectively. 


III. NONLINEAR COMPRESSION OF 
QUANTUM GREY SOLITON 


Single component 

Unlike adiabatic soliton compression, the present sys¬ 
tem takes an advantage of the exact solution to the non¬ 
linear Schrodinger equation for chirped grey soliton evo¬ 
lution. In order to see the nonlinear compression, as well 
as for experimental realization, we present the follow¬ 
ing example. For simplicity, we consider the loss/gain 
term G = 0. The Feshbach managed scattering length 
can be tuned as: as{t) = [35]. In the presence 

of an expulsive parabolic oscillator with M{t) = — A^, 
we find the chirping coefficient as: c(t) = A. Consecu¬ 
tively, all other parameters are obtained: A{t) = 
and l{t) = These solitons are necessarily chirped 

and can even be accelerated, compressed or brought to 
rest. Fig. &) depicts the nonlinear compression of the 
dark soliton; one notices that with the increasing value 
of the scattering length, the dark soliton has an increase 
in amplitude and a compression in its width. The soliton 
can be brought to rest at ^ = 0. As we increase the value 
of 0, the velocity of the soliton increases and the soliton 
propagates along the longitudinal direction. 

Therefore, the acceleration and the compression of the 
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FIG. 3. (a) The nonlinear compressin of the quantum grey 
soliton {0 — 7r/16) as it propagates. The width of the soliton 
decreases with the increasing depth, (b) The hole-like excita¬ 
tions spread out in presence of the loss/gain term (yd = 0.12). 
The other parameter values are same as in Fig. 


FIG. 4. The nonlinear compression of the dark and bright 
soliton is shown for ^ = tt/IG. As is seen, the width of the 
dark-bright solitons decreases with increasing amplitude. The 
other parameter values are same as in Fig. 0 and 0. 


soliton can be controlled through trap modulation and 
tuning of the scattering length by means of Feshbach res¬ 
onances. In the following section, we extend our analysis 
to the two component case, in order to see the effect of 
intra- and inter-atomic interactions, on grey soliton dy¬ 
namics. 


domain, where Lieb mode becomes experimentally real¬ 
izable in the background of sound mode. 

IV. LIEB-MODE EXCITATIONS 
Single component 


Binary BEC 

There are many tools, which are being used to control 
and manipulate the various parameters of solitons and 
induce changes in their shapes, which is very useful for 
developing many applications of BECs. One way is to 
vary the atomic scattering length by means of external 
magnetic field, known as Feshbach resonance. We here 
demonstrate that in the presence of an ’’expulsive” trap, 
the variation of scattering length in a particular manner 
leads to a nonlinear compression of both dark and bright 
soliton simultaneously. Fig. |^a) shows tha compression 
of the width of a dark soliton, as seen in case of single 
component. The nonlinear compression of the bright soli¬ 
ton is shown in Fig. |^b), which shows that the width of 
the bright soliton can also be compressed. As the scatter¬ 
ing length increases, the amplitude of the bright soliton 
increases. At 6> = 0, the dark soliton becomes a black 
soliton with its velocity u = 0. The amplitude of the 
bright soliton becomes maximum at this point. As we 
increase the amplitude of the bright soliton decreases 
and vanishes at 6> = 7r/2, whereas, the dark soliton takes 
a constant value. At this point, the velocity of the soli¬ 
tons achieves its maximum value, i.e., the sound velocity. 
This technique can be used to modulate the bright soli¬ 
ton into very high local density in an inverted harmonic 
oscillator potential. 

The above description for solutions of multi- 
component BEC now provides us the platform to investi¬ 
gate different modes associated with this system. In the 
following section we investigate the dispersion relations 
of Lieb and Bogoliubov modes carefully and identify a 


To study the regime of Lieb-mode, one requires the 
calculation of canonical energy and momenta of the sys¬ 
tem. The energy of the grey soliton can be evaluated 
by subtracting the energy due to the contribution of the 
background term, such that E = W — VFo, where. 





dz. 


( 11 ) 


The energy expression kFo corresponds to the profile i/iq = 
YiSiS a kinematic phase in a trap of 
similar parameter values. Given the solitary wave profile 
of Eq. (§, Eq. 0 yields. 




cos 0 ) + 2Abuao( cos 0 + cos 0 (12) 


where, the exponential pre-factor shows the effect due 
to the loss/gain of atoms from the condensate. The first 
term in the expression represents the contribution from 
the solitonic energy. The term thereafter accounts for 
the oscillator and chirp contributions; it is singular at 
0 = 7r/2. The third term, linear in 6(t), arises from 
the coupling of BEC momentum with the phase of the 
soliton.The last term being quadratic in 6(t), represents 
BEC translational energy. One must recall that we have 
started with a very general structure, keeping several 
equation parameters as time dependent. If one looks at 
the Riccati equation and considers c(t) as independent 
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FIG. 5. Dispersion relations for the soliton (blue dashed line) 
and Bogoliubov mode (pink dotted line), (a) The behavior 
of Lieb and Bogoliubov modes are shown at t = 0. Inset 
reveals that the periodicity differs when t 7 ^ 0. (b) The Lieb 
mode dispersion for 6 = 0, m presence of harmonic trap. Inset 
depicts the mode behavior, when 6/0. The other parameter 
values are same as in Fig. 0 - 


of time, then we obtain a condition: + M = 0. This 


removes the energy divergence in Eq.(12) and it exactly 


matches with [la US]. This situation we have already 
noted in our preceding section and showed that the sub¬ 
sequent expressions of the equation parameters are in 
accordance with literature [35]. Thus, through a con¬ 
trolled modulation of the external magnetic field, we are 
actually able to take into account the external trap and 
extract the usual result through a complete analytical 
model. However, the temporal dependence of the chirp¬ 
ing also carries some significant information, which we 
will elaborate after presenting the canonical momentum 
calculation. 

The canonical momentum of the solitary wave is given 

by, 

P = -i [ iP*^dz 


dz 


= e 


CsCcfq tt^ — sm20 — 20]— 2b(aocos0(13) 
u 


The first term can be attributed to soliton, while the 
second term arises from the chirped momentum. The 
energy and momentum are, respectively, normalized by 
Eq = ^0 = ^oCsCToC- is worth mention¬ 

ing that the soliton velocity can also be computed from 
the hydrodynamic relation: ^ = Au^ which matches 
with the earlier obtained results. In the limiting case, i.e., 
when the trap is switched off, all the expressions match 
with the known results nsms]. Eq. ( [T^ implies that 
a maximum momentum of Pmax = {CsCcfo^ — 2hC,(Jo) 
is obtained ioi u = 0. The Lieb mode terminates at 
P = Pmax- The Bogoliubov or sound mode is calculated 
using the usual prescription, by the application of a per¬ 
turbation and subsequent linearization of the perturbed 
Gross-Pitaevskii equation. 

We present different consequences of dispersion in 
Fig. ([ 5 ]). In Fig. ([5^) the dispersion associated with the 
grey soliton profile with well known 27r periodicity with 
t = 0 and without the trap. However, when t / 0 in 


the absence of trap, the characteristic periodicity of the 
solitonic energy will change. We have shown this behav¬ 
ior in the inset of Eig.([^). Now we would like to draw 
the attention to the Eig.([^). In Eig.([^) the solitonic 
energy shows a divergence. One may try to avoid this 
situation in experiment, we show here that this can offer 
new possibilities as well. We can clearly see that in a 
small momentum window, the solitonic energy is much 
higher than the sound mode, implying that it is possi¬ 
ble to observe this solitonic mode in a restricted region 
with controlled chirping (which can be controlled exper¬ 
imentally) before the system gets destabilized. Eurther, 
the presence of the translational motion of BEC consid¬ 
erably affects the grey soliton and one can observe that 
in a certain low momenta region, the soliton’s momen¬ 
tum can become negative, as depicted in the inset of Eig. 
#), thereby suggesting a region of negative soliton mass. 
More precisely, taking into account the momentum con¬ 
servation, the time-dependent momentum term in chirp¬ 
ing b{t) is expected to get strongly correlated with the 
COM of the soliton. This strong coupling between the 
chirped momentum and COM enables us to achieve the 
regime of negative mass. The strong coupling between 
the chirp momentum and COM can be understood from 
the effective mass profile. Therefore, we obtain the effec¬ 
tive mass, as shown in [4811^ . 


rrieff = P/V = moE{0,u), (14) 


where, ttiq = HocroC is a scaling factor. The factor E(^, u) 
can easily be obtained using Eq. (13). The effective 
mass profile is shown in Eig. Eor positive cur¬ 
rent conservation yields, to be negative for hole¬ 

like excitations. This accumulates a net phase difference 
(negative) as the grey soliton starts propagating. The 
combine effects of this negative phase accumulation and 
chirped momentum enable us to enter into the negative 
mass regime. The variation of the effective mass as a 
function of 0, as shown in Eig. (|^) provides a better un¬ 
derstanding. It shows that small values of 6>, resembles 
the small momentum values, the effective mass becomes 
negative. This is akin to the scenario depicted in the 
inset of Eig. ^jp). Eig. ^jp) replicates the same as a 
function of the inverse velocity We observe that 

the effective mass becomes negative for small values of 
u and the slope (negative) gives rise to the momentum. 
Hence, it is possible to achieve the negative mass regime, 
which emerges in a small momentum window. We em¬ 
phasize that this regime is possible to achieve experimen¬ 
tally though the chirp management within the small mo¬ 
mentum window before the system gets destabilized. In 
this context, it is worth mentioning that negative mass 
regime for gap solitons has already been proposed [50| 
and observed experimentally m- 


As a next step we show that similar dispersion mecha¬ 
nism also arises in the two component case, which reveals 
much richer physics, due to the presence of the inter¬ 
species interaction. 
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FIG. 6. The effective mass profile is shown as a function of 
(a) Mach angle 0 and (b) velocity u. This is akin to the 
scenario of negative mass in the low momenta regimes, (b) 
shows the emergence of negative mass regime for small values 
of u and the negative slope gives rise to the momentum. The 
parameter values are same as in Fig. 0- 


Formation of bound state 


The energy of the grey soliton in case of TBECs can be 
obtained from Eq. (11), where, the background energy 
Wq is obtained by assuming a constant solution. Unlike 
the single component case, we consider here b{t) = 0. 
Since, we are interested in the Lieb-mode excitation, we 
concentrate only to the energy of the quantum grey soli¬ 
ton. The parameter values, used for the constant solution 
are same as the previous case. Using the solitary (grey) 
density profiles in Eq. 0 and ( [Tq| ), Eq. ( pT| yields. 




cosO \ -KaO-^A^ cos^ 0 


-^Ci^abCrob^A^{t) cosO. 


(15) 


Similar to the single component case, the contributions 
from the harmonic trap and chirping to the energy turns 
out to be singular at 7r/2. The third and fourth terms 
are due to the presence of intra- and inter-species inerac- 
tions, respectively. The canonical momentum of the grey 
soliton is found to be. 


u 

Pa = C^oA(t)us(7r -—r — 20 — sin26>) + 2(aoc(t)l(t) cosO. 

\u\ 

(16) 

The first term in the momentum expression is the mo¬ 
mentum of the soliton and the second term is due to the 
presence of chirping. The effect of the inter-species inter¬ 
action on dispersion profile is investigated and is shown 
in Eig. 0 . At low momenta, the energy is found to be 
very high. The fact that the presence of trap consider¬ 
ably affects the dynamics of grey soliton and does not 


allow a uniform density, which leads to the strong inter¬ 
action between these hole-like excitations and trap geom¬ 
etry. The Bogoluibov (BG) mode is shown as red dots. 
In the focusing case (attractive interaction), the energy 
of these hole-like excitations lies above the single com¬ 
ponent case, which shows degeneracy (the filled squares) 
{f^ab = —0.5). As we decrease the interactions strength 
{i^ab = —0.2), the energy gets reduced, but maintains 
the degeneracy (the filled diamonds). At K.ab = 0, we 
retrieve the dispersion behavior of the single component 
case (the filled upward triangle). As we move towards 
the defocusing case of a positive sign (repulsive interac¬ 
tion), the inter-species interaction {bZab = 0.2) leads to 
a reduction of energy at each point of high momentum, 
shown as inverted triangle. If the values of the Kab is in¬ 
creased further, energy becomes negative, that indicates 
the formation of a bound states, as depicted in the figure 
by empty circles. 



• BG mode 
■ ^ab= -0-4 

*^ab= -0-08 

^1D BEG 

▼ ^ab“ 0.08 

n /Cab= 0.4 


FIG. 7. The figure shows the dispersion relation of the grey 
soliton for different values of inter-species interactions. As 
shown, for attractive (repulsive) interaction the dispersion lies 
above (below) the dispersion of the single component one. 
The other parameter values are same as in Fig. 0 and 0 - 


The intriguing aspect of bound state formation be¬ 
tween the binary components as manifested in Fig. 
can be explained in the realm of low energy scattering 
theory [52]. Here, the effective interaction between the 
two species, hZab actually relates the transition of the sys¬ 
tem from weak coupling to the strong coupling regime. 
In this transition the attractive inter-species interaction 
hiab points to the weak coupling. The gradual increase in 
the bare coupling strength results in the effective interac¬ 
tion to be positive, thereby indicating a strong coupling 
regime which, as in other cases (e.g composite boson for¬ 
mation), favors the bound state formation. This nature is 
explicated in the figure, as the solitonic energy dips below 
the zero level. In this context, we would like to point out 
that such controlled bound state formation may be useful 
for information storage and retrieval. It is worth noting 
that grey soliton excitations have non-trivial phase differ¬ 
ence; the phase is intensity dependent, whereas in case of 
bright soliton, it is of kinematic origin. The relative phase 


















can code information, which may be stored in the static 
configuration, either in the bound state or in the plane 
static configuration. The phase can be retrieved later by 
separating them and making them dynamical. Informa¬ 
tion storage and retrieval in nonlinear media, particularly 
in active atomic media has already been demonstrated in 

m- 

V. CONCLUSION 

In this article, we study inhomogeneous multi compo¬ 
nent BEC with special attention to the dispersion mech¬ 
anism. We present an analytical model to tackle the in¬ 
homogeneity and relating them to he homogeneous case. 
Later, we show that our model is general enough to repli¬ 
cate well known dispersion diagram while we model the 
interaction in a certain way [35]. The specific control 
on the external magnetic field (effectively controlling the 
interaction) automatically modifies the chirping. Inter¬ 
estingly a generic chirping (with which we started) can 
enable us to view the hole-like excitation mode exper¬ 
imentally in a small momentum window before getting 
destroyed. Our study affirms the recent description of 
the finite lifetime soliton in a trap [53] in a situation when 
we employ coherent control on the chirping. The mod¬ 
ulated chirping enforces instability for the soliton but it 


also provides a scope to experimentally observe the Lieb 
mode in a tiny momentum window. Further, we note 
that, through our formalism, it is also possible to get 
back the usual dispersion diagrams akin to the general 
understanding by means of chirp-engineering. This opens 
up the possibility for the coherent control and manipula¬ 
tion of these hole-like excitations at low energy, through 
the temporal modulation of the scattering length, as well 
as the trap. Additionally, we have shown a wide class 
of exact solutions for the coupled quasi-one dimensional 
GP equations in a trap, describing the dynamics of bi¬ 
nary condensates. We concentrated on the case where 
one component possesses a hole-like excitations and the 
other being a bright soliton. The formation of bound 
state in repulsive inter-species interaction is shown and 
discussed that it is manifested from the low energy scat¬ 
tering theory. 
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